Application of a parallel hybrid iV— body SCF code 
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Q\ ' 1. Introduction 
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The stellar dynamics of galactic nuclei are in the focus of astronomical interest. 
Q \ Observations show that many if not all galactic centres harbor supermassive 

^ | ■ black holes. Their mass determination, however, relies on stationary modelling, 

which does not provide a self-consistent picture of the evolution of galactic nuclei. 

Another point of interest is the behavior of binary black holes within a 
dense central stellar cluster. To what extend does star gas or dynamical friction 
J> ' influence the results? These applications need a high precision collisional code, 

which can handle a great number of particles (Quinlan & Hernquist 1997). 

The iV 2 dependency of computational time restricts high precision direct 
iV-body calculations to a comparatively small number of particles. In contrast, 
codes which calculate the mean potential of particle distribution do not have this 
problem but are not sufficient to treat the dynamics of central stellar clusters. 
Since many stellar systems show a clear core halo structure, one can profit from 
the advantages of both approaches in a hybrid code. The concepts of such a 
hybrid scheme and its results are described here. 

o 

H 

to ! 2. Merging of NBODY6++ with SCF 

The primary iV-body integrator for collisional systems like globular clusters is an 
Aarseth NBODYx code (Aarseth 1993). The most complex one, NB0DY6++, 
is highly efficient on serial computers. This code has been ported to parallel 
■ computers (Spurzem & Baumgart 1997) where it can show its advantages com- 

pared to iV 2 codes running on special purpose computers. Nearly all parts of 
the code, like the ones for binary or chain treatment, can run in parallel as long 
as the problem size is large enough. 

To embed NBODY6++ in a living, nearly spherical halo, we choose the 
SCF method by Clutton-Brock (1973) and Hernquist & Ostriker (1992). The 
parametrised a-law basis functions published by Zhao (1996), which are used 
in this code instead, provide a much higher flexibility. The Hermite integration 
scheme in NBODY6++ needs the input of the forces aj and force derivatives a;. 
Unlike in other mean field methods, these can be added to SCF analytically. 

The angular momentum of the particles is the criterion to distribute the 
work over both codes. Particles with a low angular momentum are more reactive 
to their neighbors and therefore given to NBODY6++. In a Plummer potential, 
this criterion assigns all spatially central particles to direct integrator. 

The SCF method relies on an analytically given distribution function for 
density and for potential. Thus, a cut off potential is very hard to handle. This 
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problem can be solved by calculating the the forces of the whole particle set 
using SCF. NBODY6++ and SCF re-calculate the forces in the subset below the 
critical angular momentum. In this process the direct code takes the difference 
force of the two field method runs as an external potential. The time steps of 
the field method are usually longer than the ones of NBODY6++. 

3. Testing 

One important test of the code is the benchmarking. The table below shows 
the computation time in seconds for 40 SCF time steps for various numbers of 
processing elements. About half the number of particles were below the critical 
angular momentum. 
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One realization of a Plummer model was evolved by the hybrid code and a 
HARP version of NBODY6++ for several crossing times in order to check the 
difference in relaxation effects. In the field method part of the hybrid code one 
observes less relaxation than what the direct method finds. Introducing Spitzer 
relaxation to SCF, however, might be a way to tackle this problem. 
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Relaxation of aspherical galaxies 
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Abstract. Colour (or, metalicity) gradients are commonly observed in 
elliptical galaxies (eg, Danziger et al. 1993). These may be attributed 
to stellar evolution processes and gas recycling, or the conditions at the 
time the galaxy formed. A tractable approach to galaxy formation is dis- 
sipationless cold collapse. This violent but incomplete relaxation (eg van 
Albada 1982) means that traces of the initial orbit distribution remain, 
and so gradients in galactic properties. Here we seek out a possible re- 
lation between orbit survival and the on-set of radial orbit instabilities 
(ROIs). With this in mind, we apply Newtonian gravity to collapsing cold 
spheroids of point masses as a model for the Lin-Mestel-Shu (1965) solu- 
tion for dust. Below we setup a test problem which we solve numerically 
and discuss. 



1. N-body integration 

Our toy model consists of a uniform sphere containing two mass components, 
one inside, the other outside, r = tt; the bounding radius = 5. Each mass bin 
is sampled with 5,000 bodies. Making mi/ni2 = 1/3, the mean density is ev- 
erywhere the same and hence the radial free-fall time is unique. With the more 
massive stars confined to it < r < 5, the intention is to favour mass ejection 
from, or otherwise disturb, the lighter central distribution. This is a robust test 
of orbit survival in energy space under violent relaxation. 

The 10,000-particle simulation was performed with the code NBODY1 on 
HARP computer. The length I of the smoothed potential = 3 x 10 -4 (in system 
units), well below close-encounter impact-parameter values. The simulation ran 
for a total of 6 initial crossing time (= t cr ). A small velocity was given to 
each star such that the virial ratio Q = —T/W = 1/200 initially. Two more 
simulations were performed with the same random seed distribution, but now 
flattening the z-axis by factors of 2 (hamburger) and 10 (pancake): this covers 
an appropriate range of spheroid aspect ratios. 

The two-body relaxation timescale for this system is ~ 60i cr (Chernoff &: 
Weinberg 1990; Spitzer 1987) so evolution over 6i cr is largely collisionless. 



1 Now at: Astronomisches Rechen-Institut, Monchhofstrasse 12-14, Heidelberg D-69120, 
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2. Discussion 



We measured orbit-mixing by monitoring in time the ranking TZt,i of each (i th ) 
star relative to one another, introducing 



Completely mixed orbits give A ~ 57.7 when S the number of rank bins = 200. 
Both spherical and 'pancake' collapses erase memory of the initial vertical (z) 
structure rather well, however this is not so for the 'hamburger' simulation: 
we find post-collapse A z s=s 42, well below the expected value for fully-mixed 
systems. Inspection of the principal components of the inertia tensor showed 
that both mass components relaxed to a triaxial figure, with axis ratios 1.00 : 
0.87 : 0.80, in rough agreement with values derived from a scaling argument 
applied to LMS-type collapses (see Boily et al. 1998). The heavy masses filled a 
volume twice as large (larger I xx per unit mass) so the lighter population orbited 
in the background potential of the heavies. 

For the case of the spherical collapse, we found at the bounce and there- 
after that the lighter masses formed a spheroidal mass profile, while the heavier 
particles remained spherically distributed. Furthermore, they were also (mildly) 
more centrally concentrated. As we let this system evolve up to 50t cr , we found 
that orbit-orbit interactions with the heavy population reestablished sphericity 
(see also Theis Sz Spurzem 1998): in the process heavy and light stars moved out- 
and in-wards respectively, opposite the trend expected of two-body relaxation. 
Near the end of the run, both populations filled the same spherical volume. 

These results suggest further investigations into processes driving ROIs (eg, 
Aguilar & Merritt 1990) to establish long-term effects of a bias in the initial mass 
distribution. 
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Simulations of star clusters and Galactic Nuclei: rotation 
and other challenges 
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1. Introduction 



With strong efforts - both observationally and theoretically - it has now become 
possible not only to find evidence for the existence of black holes in the centres 
■^j- ' of galaxies, but also to set up an expression for the mass distribution of black 

holes with respect to the mass of the parent galaxy - or its spheroid component 
(Kormendy & Richstone 1995). In deriving quantities of central objects from 
the orbit structure of the surrounding star clusters it is crucial whether 2- or 
qq \ 3-integral models are made and which integrals are chosen. Using E and J z van 

■ der Marel et al. (1998) analyze the central cluster of M32 and find significant 

amount of rotation around a black hole mass of roughly = 3.4 x 1O 6 M0. 
The relaxation time at least in this system is several times less than its age, so 
\ that corresponding processes may play an important role in the evolution and 

build-up of central dark objects and their surrounding clusters. 

Only recently first efforts are being started to model rotating galactic nuclei 
including the effects of stellar two-body relaxation, star accretion onto the black 
hole, movement and growth of one or two central black holes. There is an 
iV-body approach by Arabadjis & Richstone (1998ab) using a TREE-code. We 
have undertaken an effort to use a simplified (2 integrals of motion, equal masses) 
2D-Fokker-Planck model to follow the evolution of self-gravitating, rotating star 
clusters into core collapse without central mass (Einsel & Spurzem 1998). The 
/S | gravo-gyro catastrophe, earlier postulated by Hachisu (1979) by using a linear 

stability analysis, has been rediscovered, followed by a self-similar contraction 
phase. 

Now, as a second step, we study multi-mass star clusters. The well-known 
mass segregation instability is now accompanied by the gravo-gyro and subse- 
quent effects governed by the interplay between rotation and two-body relax- 
ation. 

In the Figure the evolution of the central angular velocity of the different 
stellar masses (indicated at the lines) is plotted as a function of total central den- 
sity (which can be considered as a scaled time). The initial model was a Michie- 
King like model with Wo = 3, dimensionless rotational parameter ujq = 0.30 and 
a Salpeter mass spectrum. We observe that a fast rotating core of the heavy 
masses is built up, and conclude that the hitherto not carefully analysed effect 
of rotation in relaxing star clusters introduces unexpected and very interesting 
kinematic features. 
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Figure 1. Evolution of central angular velocities versus centralden- 
sity of all 10 mass bins used in the calculation for an initial model with 
Wo = 3 and ujq = 0.90. The highest mass bin with m = 1.5 Mq is the 
uppermost curve, while the lowest mass bin with m = 0.15 Mq can be 
found as the lowest curve. 



This models will be published in more details elsewhere (Einsel & Spurzem 
1999), and the inclusion of a central massive star-accreting black hole is in 
planned. 
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